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We propose an exactly solvable model to reveal the physics of the interplay between interaction 
and disorder in bosonic systems. Considering interacting bosons in a double-well potential, in which 
disorder is mimicked by taking the energy level mismatch between the two wells to be randomly 
distributed, we find "two negatives make a positive" effect. While disorder or interaction by itself 
suppresses the phase coherence between the two wells, both together enhance the phase coherence. 
This model captures several striking features of the disordered Bose-Hubbard model found in recent 
numerical simulations. Results at finite temperatures may help explain why a recent experiment 
did not find any evidence for the enhancement of phase coherence in a disordered bosonic system. 



Physics of disorder in both bosonic and fermionic 
systems has been attracting great interest from physi- 
cists for decades. While the effect of disorder on non- 
interacting particles can be well described by Anderson 
localization^ , understanding the interplay between in- 
teraction and disorder still remains challenging 2J . Re- 
cent developments in studying cold atoms provide new 
opportunities to reveal the nature of interacting quantum 
particles in the presence of disorder j3]. In atomic sys- 
tems, the strength of both interaction and disorder can 
be well controlled experimentally, unlike in solid state 
systems. There has been significant recent progress in 
studying disorder effects in both non-interacting and in- 
teracting atomic systems [4HT0]. 

From the theoretical side, to understand even the spin- 
less bosons in the presence of both interaction and disor- 
der is a non-trivial problem. Scaling analysis, renormal- 
ization group theory and sophisticated numerical simula- 
tions are often used |llffi"8] . It is desirable in this context 
to have some exactly solvable models, which capture all 
the ingredients, such as interaction, disorder and finite 
temperature effects. In this letter, we propose and solve 
such a model, which makes it possible to reveal the un- 
derlying physics transparently. Our minimal model cap- 
tures the interplay between interaction and disorder in 
bosonic systems, showing conclusively that weak disor- 
der in the presence of interaction generically enhances 
the phase coherence of the system. 

The model we propose describes interacting bosons in 
a disordered double-well potential. The Hamiltonian can 
be written as 

H = -t{b\b R +c.c)+-^n a {n a -l)+^{n R -n L ), (1) 

a 

where b L (br,) and b R {b R ) are the creation(annihilation) 
operators in the left and right well respectively, a — L,R, 
n L = ^iphinR = b R bn) is the number operator in the left 
(right) well, t is the tunneling amplitude between the two 
wells, U > is the onsite interaction, and e is the energy 
level mismatch between the two wells. A special case of 



this model (e = 0) has been studied both theoretically 
and in cold atom experiments 20 22]- In our case, e is 
randomly distributed according to a certain probability 
function P(e), thus simulating disorder. In experiments, 
this can be realized by randomly tilting the double-well 
potential. The thermodynamic quantities will be the en- 
semble averaged values. In this paper, we will focus on 
the case P(e) = 1/(2A), e € [—A, A], where A character- 
izes the disorder strength. Other distribution functions 
for e do not change the qualitative conclusion presented 
in this paper. Eq.Q can be viewed as a two site version 
of the extensively studied disordered Bose-Hubbard(BH) 
model with random onsite energies. 

Despite the simplicity of this model, many interesting 
phenomena are found by solving this model with N par- 
ticles. An illuminating example is the case e = 0. It 
has been shown that the ground state has a crossover 
from the coherent state to the Fock state when U/t 
increases [22] . which mimics the well-known phase transi- 
tion in the thermodynamic limit from a superfluid con- 
densate to a Mott insulator state in the full BH model. 
The physics for either the crossover or the phase tran- 
sition between the phases is the same in the minimal 
model and the BH model, namely, interaction suppresses 
number fluctuations and spatial phase coherence. In the 
case where e is randomly distributed, we will see that 
the minimal model qualitatively captures the physics of 
interacting bosons in a disordered potential for a number 
of fundamental questions, even though it can not answer 
some quantitative questions, such as where there is a di- 
rect transition from Mott insulator to superfluid. 

Before we discuss the details, we briefly summarize 
the main qualitative questions that we are going to an- 
swer in this paper. Ql: Whether disorder can enhance 
the phase coherence between the two wells, defined by 
C = {b\bpi)/N = {b\bL)/N, in some parameter regimes, 
in contrast to the intuition that disorder must always de- 
stroy the phase coherence? Q2: Whether interaction can 
also enhance the phase coherence in the presence of dis- 
order, though without disorder it is known that repulsive 
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interaction usually suppresses phase coherence? Q3: In 
the clean system, e = 0, the compressibility of the sys- 
tem vanishes when the phase coherence is destroyed by 
interaction. Is this still true in the disordered case? Q4: 
How does the phase coherence between the two wells de- 
pend on the temperature? We will see that the answers 
to these four questions provide insight into some striking 
features of the quantum phase diagram of the disordered 
BH model. Our answers to these questions on the min- 
imal model also shed light on a recent disagreement be- 
tween the prediction of numerical simulations |17U18j and 
experimental observations [5]. 

The answer to Ql: To solve the disordered prob- 
lem, we start from the case of N particles in the double- 
well with a fixed e. We write the Schrodinger equation 
H\V) = in the Fock space. Define \l) = \N L ,N R ) = 
+ where I = 0, ±1, ±2, ... ± N/2. For sim- 
plicity, we assume N is a large even number. Expanding 
|*) = Sj^lO) we obtain 

(E - Ei)\l) = -tM u+1 \l + 1) - tM t _ ltl \l - 1), (2) 

where E l = Ul 2 + d + f (A 2 /2 - N), M Ll +i = M l+U = 
v / N/2(N/2 + 1) - 1(1 + 1). The eigenenergies and eigen- 
functions can be easily computed by exact diagonaliza- 
tion. At zero temperature, the phase coherence between 
the two wells can be characterized as C e = ^(6^6/?) € = 
Jr -^M+i^VJ+i! where ^ is the ground state wave 
function in the Fock space, (O) is the expectation value 
of the operator O in the ground state, and the subscript 
implies a fixed e. The results for the disordered case can 
then be obtained by averaging those results for fixed e 
according to the distribution function P(e), 

<0) = / deP(e)(0) e . (3) 
J-A 

The results for C e for both noninteracting case and fi- 
nite U are shown in Fig.(lA). The numerical results with 
the total particle number N — 100 have been chosen to 
illustrate the physics, which does not depend on the ex- 
act value of the total particle number. When U — 0, C e 
quickly decreases with increasing e, since noninteracting 
particles can only tunnel from one well to the other if the 
energy mismatch e is smaller than the tunneling t. The 
behavior of C e as a function of e depends explicitly on U, 
as shown in Fig.l. The dependence of C t on e becomes 
more extended with increasing U, though C e (e = 0) de- 
creases with increasing U. More importantly, the curves 
become non-monotonic when the discretization of E\ in- 
duced by interaction becomes significant for U > tN. 
Local maxima of C e emerge at certain values of e. For 
those values e* = ?7(1 — 2/*), where I* is an integer num- 
ber, Ei* — Ei*_\ is satisfied. In other words, the state 
\N/2 — I*, N/2 + I*) has the same energy with the one 
\N/2 - I* + 1, N/2 + I* - 1). There is no extra inter- 
action energy cost for one particle tunneling from one 
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FIG. 1: A: C t as a function of e/tN at different interaction. 
B: C as a function of A/tN. Dashed(purple), dash-double- 
dotted(green) , solid(brown), dotted(blue), dash-dotted(red) 
represent U/tN = 0, 0.3, 1, 2.5, 9.95 respectively. Insets are 
results for C t and C for U/tN — with the same labels as 
the main figures. 

well to the other. The tunneling will be enhanced at e*, 
forming gently varied bumps on the curves for C £ . This 
non-monotonic behavior of C c becomes even more dra- 
matic when U increases. For large values of U , only in 
narrow regions near e, tunneling is strongly enhanced. 
As a result, resonance features emerge with sharp peaks 
located at e* sitting on the slowly decaying envelope of 
C e . In the large U limit, E\ > E t , = if I ^ I* 

or I* — 1, a two level approximation can be made, and 
it is straightforward to show that C e » ~ M;._i ! ;*. C e » 
decreases slowly as e* increases, as shown in Fig.(lA). 

Now we turn to the quantity C = J A A deC e , which 
is related to the area between e 6 [—A, A] below the 
curve of C € . When U is small, averaging C e according 
to Eq.([3| smooths out the small bumps on the curves. 
The non-monotonic behavior of C is largely suppressed, 
as shown by the solid(brown) curve in Fig. (IB). On the 
other hand, we have seen that the resonance feature of 
the curves becomes significant for large U. When A is 
not large enough, the contribution from regions under 
the peaks to the total area is dominant. As a result, 
peaks can still be clearly resolved for small values of A, as 
shown by the dotted(blue) and dash-dotted (red) curves 
in Fig. (IB). For large A, the peaks become wiggles on the 
top of the slowly decaying curves of C. It is also quite 
clear that C first increases when the disorder strength A 
grows from zero, if A < A* , where A* = U is the location 
of the first peak. This fact shows the significant difference 
of the disorder effect between the non-interacting and 
interacting systems. In the presence of interaction, the 
disorder effect is non-monotonic. 

It is helpful to look at the contours of C as a function 
of A/t and U/t in order to obtain a complete picture. 
Fig. (2) shows that for large U, the contours first bend to 
the right hand side when A increases from zero. We note 
that the contours of C have wiggles at large U. These 
wiggles arise from the structures in C as discussed in the 
last paragraph. For a fixed value of U, away from the 
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FIG. 2: A contour plot of C as a function of A/tN and C//tJV 
at zero temperature with the values of contours marked down. 



positions of the peaks, C changes slowly, corresponding 
to the parts of the contour which are nearly parallel to 
the A axis. When approaching the peaks, C quickly in- 
creases. As a result, the contours bend towards the U 
axis, forming a wiggly shape. The topology of Fig. (2) 
is very similar to the nontrivial structure of the phase 
boundary obtained by recent numerical simulations for 
the disordered BH model at an integer filling [T71 ITS], In 
the latter case, the phase boundary can be viewed as the 
contour for the order parameters (bi) =0 in the thermo- 
dynamic limit. 

We emphasize that the similarity of the topology be- 
tween Fig. (2) and the phase diagram of the disordered 
BH model is not accidental. The lattice model can be 
viewed as the thermodynamic limit of the two-site prob- 
lem. Moreover, the physics of the interplay between dis- 
order and interaction is the same in both cases. Inter- 
action discretizes the energies of the Fock states in each 
site, and suppresses the tunneling between different sites. 
However, disorder introduces relative energy shifts of the 
Fock states at different sites. In a randomly distributed 
disorder potential of large enough strength, there are al- 
ways possibilities for neighboring sites to have nearly de- 
generate Fock states. Effectively, the tunneling of the 
particles as well as the spatial phase coherence will then 
be enhanced. If the disorder strength increases further, 
the weight of those configurations favoring tunneling in 
all the configurations of the random potential decreases. 
The phase coherence is thus eventually suppressed by 
very strong disorder, as seen also in the numerical sim- 
ulations of the full disordered BH model [T5HT8] . It is 
worthwhile to point out that a quantitative understand- 
ing of the establishment of a long range order in the dis- 
ordered BH mode requires taking into account long-range 
correlations beyond neighboring sites[19]. Nevertheless, 
our two-site model qualitatively reveals the underlying 
physics for the enhancement of phase coherence through 
weak disorder. 

The answer to Q2: Having answered Ql, the an- 
swer to Q2 becomes clear. Fig. (3A) shows the depen- 
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FIG. 3: Left: C as a function of U/tN at A/tN = 0(dashed, 
purple), 20.4(dash-double-dotted, green), 55.5(solid, brown), 
83.7(dotted, blue) and 140(dash-dotted, red). Right: com- 
pressibility(solid red) and C (dash-dotted blue) as a function 
of A/tN at U/tN = 0.5. The dashed(magenta) line repre- 
sents the asymptotic value of C at large interaction limit. 



dence of C on U/t at fixed A/t. For A = 0, interaction 
suppresses phase coherence monotonically. For A ^ 0, 
interaction first enhances phase coherence before sup- 
pressing it. These results are reminiscent of a similar 
behavior of the superfluid density as a function of U/t 
in the full disordered BH model jT51H7] . In both cases, 
interaction screens the disordered potential, since some 
particles occupy the sites with lower on-site energies and 
thus smooth out the effective potential for the remaining 
particles. The spatial phase coherence is then enhanced. 
However, if interaction becomes very strong, interaction 
itself eventually destroys phase coherence or superfluid. 

We have so far seen an interesting "two negatives make 
a positive" effect on the spatial phase coherence of a sys- 
tem owing to the interplay between interaction and disor- 
der. With only disorder or interaction, the spatial phase 
coherence is suppressed either by the single particle lo- 
calization or the emergence of Mott state. When both 
interaction and disorder are present, our exactly solvable 
model clearly shows that the spatial phase coherence is 
enhanced in the parameter regime where they are com- 
parable in strength. 

The answer to Q3: To extract the compressibil- 
ity K e , we first calculate the chemical potential, /i £ = 
{E°(N + SN) - E°(N))/6N, where E°(N) is the ground 
state energy of N particles at fixed e. From the depen- 
dence of fj, e on the particle number N, n t — dN/d/i^ can 
then be obtained. The ensemble averaged compressibil- 
ity k can be calculated in the same manner as C, i.e., 
k = J_ A deK e . The dependence of k on A is shown in 
Fig.(3B). For comparison, the result of C for the same U 
is also shown. 

An interesting feature of the compressibility in the 
minimal model is that it remains constant for large dis- 
order strength, even though the phase coherence C de- 
creases for very large disorder. In the large e limit, 
almost all the particles fall into one of the two wells, 
E° -> UN(N - l)/2 - eN. It is easy to see k £ -> U~ x . 
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FIG. 4: Left: C/3 as a function of A/tN at different tempera- 
tures for a fixed interaction U /tN = 2. From top to bottom, 
T/tN = 0.4, 1, 2, 3. Right: a contour plot of C as a function 
of U/tN and A/tN at the temperature T/tN = 5. 



The ensemble averaged value k for large A is mainly de- 
termined by the contribution from n e in the large e re- 
gion. As a result, k also approaches a constant value 
U^ 1 when A is very large. We can view the state in 
this regime as an analog of the bose glass phase in the 
thermodynamic limit, which has vanishing order param- 
eter = 0(or supefiuid density p s = 0) but a finite 
compressibility n. To fully understand the bose glass 
phase in the disordered BH model, one needs to con- 
sider correlations beyond nearest neighboring sites, which 
also contribute to the finite compressibility of bose glass 
phase, despite the absence of a long range order. Nev- 
ertheless, our two-site model provides a simple example 
to demonstrate why the compressibility of a disordered 
system can remain finite after the phase coherence has 
been destroyed at large disorder strength. 

The answer to Q4: We have seen that increasing dis- 
order at large interaction leads to an increase of phase co- 
herence at zero temperature. The topology of the phase 
diagram at zero temperature obtained from recent nu- 
merical simulations of the full disordered BH model also 
indicates a transition from the insulating phase to the su- 
perfluid phase with increasing disorder strength |17l I18j . 
However, a recent experiment did not find any evidence 
for the increase of the condensate fraction with increasing 
disorder strength 9J. We will show that the topology of 
the contours of C changes at finite temperatures, namely, 
the contours may bend to the left hand side when the dis- 
order strength increases from zero. This fact may help 
understand why the experiment at finite temperatures 
did not observe disorder enhanced phase coherence. 

The value of the phase coherence at finite 
temperatures can be calculated by Cp tt = 
{Y, n e- 0En M u+1 ijjfipf +1 )/Z, where the subscript j3 
denotes the thermal average, <f)f is the nth eigen- 
function with eigenenergy E n , Z = 53„e _ ^ B ", and 
Cp — J A A deC/3^. The results for Cp as a function 
of A for a fixed U are shown in Fig.(4A) at different 



temperatures. At low temperatures, the wiggles on 
the curve of Cp retain. When temperature increases, 
the wiggles are gradually suppressed. At high enough 
temperatures, the wiggles completely vanish, leading to 
a monotonic decrease of Cp . Consequently, the topology 
of the contours of Cp changes completely, as shown in 
Fig(4B). Since the minimal model appears to describe 
the full disordered model rather well qualitatively, it 
is reasonable to assume that the phase diagram for 
the lattice model at finite temperatures will change 
in a similar manner as Fig.(4B). It has also pointed 
out that the weak superfluidity in the so called finger 
region, corresponding to the regions near the wiggles in 
Fig. (2) in our case, can be easily suppressed by finite 
temperature effects [TTI |2"5] . 

As a conclusion, we have proposed and solved a simple 
double-well model which incorporates many key ingredi- 
ents of the disordered BH model. It strikingly captures a 
number of novel features of the quantum phase diagram 
of the full BH model. The "two negatives make a posi- 
tive" effect is expected to be a general feature when the 
strengths of disorder and interaction are comparable. We 
believe that our minimal model could be a simple theo- 
retical paradigm for understanding the details of the full 
disordered interacting quantum phase diagram in many 
situations. We expect our work will stimulate more the- 
oretical studies to go beyond the two-site model and take 
into account the long range efforts in a lattice for a com- 
plete understanding of the disordered BH model. 
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